Photoemission evidence for a Mott-Hubbard metal-insulator transition in VO2 
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The temperature (T) dependent metal-insulator transition (MIT) in VO2 is investigated using 
bulk sensitive hard x-ray (~ 8 keV) valence band, core level, and V 2p-3d resonant photoemission 
spectroscopy (PES). The valence band and core level spectra are compared with full-multiplet cluster 
model calculations including a coherent screening channel. Across the MIT, V 3d spectral weight 
transfer from the coherent (d 1 C final) states at Fermi level to the incoherent (d°+d 1 L_ final) states, 
corresponding to the lower Hubbard band, lead to gap-formation. The spectral shape changes in 
V Is and V 2p core levels as well as the valence band are nicely reproduced from a cluster model 
calculations, providing electronic structure parameters. Resonant-PES finds that the d x L_ states 
resonate across the V 2p-3d threshold in addition to the d° and d 1 C states. The results support a 
Mott-Hubbard transition picture for the first order MIT in VO2. 



PACS numbers: 79.60.-i, 71.30.+h 

VO2, a d 1 electron system, exhibits a sharp first-order 
metal-insulator transition (MIT) as a function of temper- 
ature (T), at T MI = 340 KX The high-T metal phase has 
a rutile (R) structure, while the low-T insulating phase 
has a monoclinic (Mi ) structure with zig-zag type pairing 
of V atoms along the c-axis.— Magnetically, the metallic R 
phase shows enhanced susceptibility (x) with an effective 
mass m*/m ~ 6, while the insulating M\ phase is non- 
magnetic. Consequently, VO2 has attracted enormous 
attention in terms of a Mott-Hubbard (MH) correlation- 
induced versus a structural Pcierls-type MIT. 

Goodcnough proposed that the lowest energy t2 g states 
split into the d\\ band and the 7r* bands. In the R phase, 
the d\\ band overlaps the tt* band.— In the Mi phase, 
the zig-zag type pairing consists of (i) a tilting of the 
V-V pairs which lifts the tt* band above the Fermi level 
(E-p), and (ii) the d|| band opens up a gap because of 
the splitting into a filled bonding and empty antibond- 
ing band caused by V-V pairing. Zylbcrstcjn and Mott 
attributed the MIT to an on-site Hubbard type Coulomb 
energy U in the dn band^ Pougct et al. showed that pure 
VO2 under uniaxial stress, and Cr-doped VO2, exhibit an 
insulating (M2) phased in which, half of the V atoms 
form pairs, and the other half form zig-zag chains that 
behave as spin- 1/2 Heisenberg chains in NMR and EPR 
experiments.— Another intermediate phase [triclinic (T r )] 
connecting M% — > Mi , showed a bonding or dimerization 
transition of the Heisenberg chain V atoms, with a con- 
comitant decrease of \- Thus, it was concluded that the 
Mi insulating phase of VO2, as well as the M2 and T r 
phases, arc MH insulators. In an alternative picture, VO2 
was explained in terms of a Peierls insulating phase due 
to the structural phase transition. The electron-phonon 
interaction was considered to be important from Raman 
scattering^ and x-ray diffraction studies^ An ab-initio 



molecular dynamics study in the local-density approxi- 
mation (LDA) described VO2 as a charge-ordered band 
(Peierls) insulator^ but a genuine gap was not obtained. 
Also, the linear resistivity of metallic VO2, analysed in 
terms of an electron-phonon effect in a Fermi liquid, gives 
an invalid mean free pathi^ The LDA approach for VO2 
was discussed in comparison with Nb02 and M0O2, but 
does not provide a gap in the density of states (DOS).— 
Combining O Is x-ray absorption spectroscopy (XAS) 
and ultraviolet (UV)-photoemission spectroscopy (PES), 
the value of the d\\ band splitting was estimated to be 
2.6 eV, whereas the LDA calculates it to be 2.0 eV. The 
authors concluded that the structural distortion is dom- 
inantly responsible for the d\\ band splitting in the insu- 
lating phased 

The gap in the DOS, and suppressed moment, are cor- 
rectly obtained in the cluster dynamical mean field theory 
(c-DMFT) calculations of Biermann et al., showing the 
crucial role of strong Coulomb interactions and structural 
distortions in V02*^ The c-DMFT gives an isotropic 
metal, with the d 1 electron occupancy: 0.36 in the d|| 
(= dig) band and 0.32 each in the tt* (= e g doublet) 
bands. In the insulating phase, occupancies change to 0.8 
in the du (= ai g ) band and ~ 0.1 in each of the 7r* (= e g 
doublet) bands. Biermann et al. conclude that VO2 ex- 
hibits a correlation assisted Peierls transition. Haverkort 
et al., using polarization dependent XAS, have shown 
the importance of orbital occupancy switching in the V 
3d states across the MIT^ 3 - in very good agreement with 
the occupancies calculated in a cluster model, as well as 
c-DMFT results. Haverkort et al. conclude an orbital 
assisted collaborative Mott-Peierls transition for VO2. 

VO2 has provided other important results: an electric- 
field-induced MITii and a photo-induced MIT^ The 
electric-field-induced MIT, in the absence of a structural 
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transition, favors a MH transition. From the ultrafast (< 
psec.) time scale of the photo-induced MIT in V0 2 r^ the 
study concluded that the accompanying structural tran- 
sition may not be thermally initiated. In V02/Ti02:Nb 
thin films, the photo-induced electron- hole pairs result 
in a reduced resistivity by 10 -3 times in the insulating 
phase and a surface photovoltage effect^ We recently 
reported the photo-induced electronic structure changes 
using soft x-ray(SX)-PES, showing the combined effects 
of rigid band behavior and correlation effects*^ Several 
PES studie a 17 ' 18 ' 19 i 20 ' 21 i 22 have addressed the electronic 
structure changes across the T-dependent MIT in VO2. 
Early work discussed the role of strong correlations versus 
structural effects in VOiM^o. A T-dependent UV-PES 
study with an analysis of surface and bulk components, 
suggested the importance of both electron-electron and 
electron-phonon interactions for the electronic structure 
of VO2 3i More recently, infrared spectroscopy and nano- 
scale microscopy was used to reveal a divergent quasi- 
particle mass of the metallic carriers an approaching the 
insulating phased 3 - The insulating phase of V0 2 was clas- 
sified as a Mott insulator with charge ordering. Thus, we 
felt it important to study the MIT in VO2 using bulk 
sensitive hard x-ray(HX)-PES, complemented by full- 
multiplct cluster calculations and V 2p-3d resonant-PES. 
Recent studies using HX-PES have shown new features in 
core level spectra which could be consistently explained 
in terms of well-screened peaksi 24 i 25 ' 26 A study on va- 
lence band and core levels of V 2 03 using HX-PES has 
shown a clear correlation between the well-screened peak 
and the coherent peak at Ep£Z- Theoretical studies have 
also shown the importance of coherent screening 2 ^ as well 
as non-local screening 2 ^ in core level spectra. 

In this work, we study VO2 thin films using HX-PES 
(hu = 7937 cV). Since the mean free path (MFP) of 
electrons using HX-PES is significantly longer (6 - 10 
nm) than that using UV or SX-PES (< 2 nm)^ it pro- 
vides the bulk electronic structure. However, the pho- 
toionization cross sections (PICS) become very small (~ 
1CP 5 times) compared to UV or SX photons and severely 
limit valence band studies. In order to enhance signal 
intensity, and since HX attenuation lengths are much 
longer than relevant MFPs, we adopted a grazing inci- 
dence geometry^ This necessitated the use of VO2 thin 
films grown on a single crystal surface of TiC^Nb, allow- 
ing us an incidence angle of about 1°. We measure the 
valence band, the V Is, V 2p and O Is core level spectra 
as a function of T. Resonant-PES across the V 2p-3d 
threshold, and a cluster model calculation of the valence 
band and the core level V Is, V 2p spectra provide a 
consistent picture of electronic parameters for V0 2 . 

A 10-nm-thick V0 2 film was epitaxially grown on the 
(001) surface of 0.05 wt% Nb-doped Ti0 2 single crystal 
substrates using the pulsed laser deposition technique, 
as described elsewhere.— The tensile strain results in a 
first order hysteretic transition just like bulk V0 2 , but 
at a slightly lower temperature of Tmi ~ 292 K. HX- 
PES and SX-PES experiments were carried out at un- 
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FIG. 1: (Color online) (a) Valence-band spectra measured 
across the MIT using hard x-ray {hv = 7937 eV), compared 
with the cluster model calculation. Inset shows an expanded 
view (0-3 eV) of the V 3d states near Ep. (b) The main 
final state configuration spectra in the metallic phase and the 
insulating phase. 



dulator beam line BL29XUL, SPring-8 (Ref. |H) using 
a Gammadata-Scienta R4000-10 KV spectrometer— and 
BL17SU, SPring-8 using a Gammadata-Scienta SES2002 
spectrometer^ The energy resolution was set to 200 ~ 
250 meV. The measurements were carried out in a vac- 
uum of 5 x 10~ 8 Pa, at 270 and 320 K. The E F of V0 2 
was referenced to that of a gold film evaporated onto the 
sample holder in all measurements. 
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FIG. 2: (Color online) Resonant-PES across the V 2p-3d threshold in (a) the metallic phase and (b) the insulating phase, (c) 
V 2p XAS spectra. 



Figure 1(a) shows the raw data of the entire valence- 
band of V0 2 measured using HX-PES at 320 and 270 K, 
normalized for area under the curve. The valence band 
structure consists of the dominantly V 3d band near Ep 
[0-2.5 eV binding energy] and a broad band usually con- 
sidered as the dominantly O 2p band between 3 to 10 eV 
in previous results i 19 ' 21 Based on PICS and a compari- 
son with SX-PES results,— the spectral shape of the 7.5 
eV feature suggest that the V 4s contribution overlaps 
the O 2p states. The energy position and width of the 
so called O 2p band change slightly across the MIT. con- 
sistent with SX-PES results^ Figure 1(a) inset shows 
an expanded view (0-3 eV) of the V 3d states near Ep, 
which provides a direct picture of the gap formation in 
VO2. The 320-K spectrum shows a peak centered at 0.3 
eV and a clear Fermi edge, indicating the metallic state, 
and another broad weak feature centered at about 1.5 
eV. These features correspond to the coherent band at 
Ep and the incoherent (lower Hubbard) band, respec- 
tively. The larger intensity in the near Ep DOS using 
HX-PES is attributed to the bulk sensitive V 3d states, 
since the PICS of V 3d is higher than O 2p at high- 
photon (-kinetic) energy. The 270-K spectrum near Ep 
shows large changes compared to the 320-K spectrum, 
with a peak centered at 1 eV and negligible intensity at 
Ep, indicative of an energy gap of about 0.2 eV. 

Moreover, in order to investigate the origin of these 
structures in detail, we carried out calculations for a 
VC-6 cluster model in Oh symmetry, as described in ear- 
lier work. 24 ' 35 For simplicity, we assume a octahedral 
symmetry throughout. Since the crystal field effect on 
peak intensities and positions in photoemission is small 



even for the (generally dominant) cubic contribution, it 
is clear that the neglected lower symmetry terms would 
not change the spectral shapes in any significant way as 
far as photoemission is concerned. A study of lower sym- 
metry is beyond the scope of this work. In addition to 
the usual charge transfer energy, A from the O 2p ligand 
band to the upper Hubbard band, the charge transfer en- 
ergy from the coherent band to the upper Hubbard band 
is defined as A*. The ground state is described by a lin- 
ear combination of following configurations: 3d 1 , 3d 2 L, 
3d 3 L 2 , 3d 2 C, 3d 3 LC, 3d 3 C 2 . The parameters were set as 
follows: A = 4.0 eV, A* = 0.2 eV, U dd = 4.5 eV, lODq = 
1.2 eV, V = 2.4 eV, R v = 0.8 eV, in addition, U dc {ls) = 
8.0 eV, U dc (2p) = 6.5 eV, R c (ls) = 0.9 eV, and R c (2p) = 
0.8 eV for the core level calculations [Fig. 3(b) and (c)]. 
Theoretical calculations for transition metal impurities 
in a solid give a U dd ~ 7.5 eV for V, which gets reduced 
by typically ~ 60 % in an oxide^ The calculations re- 
produce the experimental spectra in the metallic and the 
insulating phases. All the parameters for the calculated 
spectra for the metallic and insulating phases are equal, 
except for the hybridization between the central V atom 
3d orbitals and the coherent band, V*. V* — 0.48 eV for 
the metallic phase, and it was set to V* = eV for the 
insulating phase. This confirms the importance of the 
coherent screening channel in reproducing the metallic 
phase of VO2 . The discrete levels obtained from the cal- 
culations were broadened by a Gaussian function of full 
width at half maximum AE = 0.35 eV. 

Figure 1(b) shows the main final state configuration 
(3d 1 £7, 3d , and Sd 1 ^) spectra in the metallic and the 
insulating phase. The total spectra are also shown in 
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FIG. 3: (Color online) (a) O Is, (b) V 2p, and (c) V Is core 
level spectra measured across the MIT using hard x-ray (hv 
— 7937 eV), compared with the cluster calculations for the 
metallic phase and the insulating phase. 



Fig. 1(a). In the metallic phase, the coherent screen- 
ing state 3d 1 C is dominant near Fermi level. On the 
other hand, the locally screening 3d state has a large 
intensity at a binding energy of ~ 1 eV with weak con- 
tribution from d l L states. In insulating phase, the co- 
herent screening states vanish and the locally screened 



3d with d x L admixture states remain around 1 eV. This 
feature corresponds to the effective lower Hubbard band, 
nicely matching the incoherent band in PES spectrum. 
The main contribution of d x L state appears around 4 eV 
in the metallic and insulating phases. We have checked 
that the ground state configuration is a mixed character 
[MH+ charge transfer (CT)] state. The difference com- 
pared to an earlier study^l stems from a dominantly MH 
description in the present study compared to dominantly 
CT description in Ref. H3. Specifically, Ref. uses A = 
2.0 eV and U = 4.5 eV, while in present case, we have 
used A = 4.0 eV and U = 4.5 eV. 



The single-site DMFT (Ref.|3J) and the c-DMFT (Ref. 
12h calculations, both indicate large changes in charge 
redistribution across the MIT. The c-DMFT calculation 
predicts a weak intensity lower Hubbard band in the in- 
sulating phase at ~ 1.8 eV, with the quasi-particle origin 
feature shifted to 0.8 eV, for a U = 4 eV and J = 0.68 eV. 
However, the HX-PES shows no feature at ~ 1.8 eV, but 
a single peak at 1.0 eV that is attributed to 3d state (the 
lower Hubbard band) in the cluster model calculation re- 
sult. The c-DMFT study points out that a lower U (e.g. 
2 eV) also stabilizes a gap, provided J is also small. How- 
ever, it probably implies that the quasi-particle feature 
will move to even smaller energies than 0.8 eV. Instead, a 
c-DMFT with a larger U and J is probably more appro- 
priate for V02- As has been discussed, while single-site 
DMFT results can also give a gap with a large U (<~ 5 
eV), it also gives local moment. While additional calcu- 
lations can provide quantitatively accurate comparisons, 
both, the single site and c-DMFT would require a U > 4 
eV. Thus we have set Udd = 4.5 eV as a reasonable value 
in present calculations. 

Using XAS and resonant-PES across the V 2p-3d 
threshold, we show the enhanced V 3c? electronic struc- 
ture. It resonates due to interference between the direct 
channel (p 6 d n + hv — > p^d' 1 ^ 1 + e) and the photoabsorp- 



tion channel (jfd n + hi 
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citations. Resonant-PES is well-known for probing this 
behavior in correlated oxides like NiO, CuO, etc. and /- 
electron systemsi 39 i 40 Using photon energies marked in 
the XAS profile of Fig. 2(c), we measured valence band 
PES in the metallic [Fig. 2(a)] and the insulating phase 
[Fig. 2(b)]. These spectra are normalized to the scan 
time and the incident photon flux. The off-resonance 
spectrum is very similar to the SX-PES spectrumi 17 ' 18 i 22 
The spectra show a typical resonance enhancement of 
the 0-2 eV feature as a function of hv : with a maxi- 
mum at hv = 518.2 eV (£3 peak in XAS). This con- 
firms 3d electron character of the 0-2 eV feature. More- 
over, the higher-binding feature of O 2p band gets en- 
hanced at a lower energy (~ 514 eV), and is further 
enhanced at 516 eV. At higher hv, a strong intensity 
Auger feature shows up at higher binding energy and 
tracks the increase in hv [marked in Fig. 2(a) and 2(b)]. 
The valence band calculations suggest the existence of 
the main d x L state overlapping the O 2p band. It res- 
onates due to interference between the direct channel 
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(d 2 L + hv — > d x L + e) and the photo-absorption chan- 
nel (p 6 d 2 L + hv -> p 5 d 3 L -> p^L + e). Note that the 
c^L final state position is not clearly fixed in the present 
calculation, because the present cluster calculations give 
a discrete level for the d 1 ^, unlike actual d^L states hav- 
ing a band like energy distribution. The Auger features 
are hardly observed in the 0-2 eV resonant feature, very 
similar to the case of Ti2 03 resonant-PES, another 3d 1 
system4i 

Figure 3(a) shows the Ols core level HX-PES spectra 
at 320 and 270 K. The O Is core level shows clear spectral 
shape change; a symmetric peak in the insulating phase 
transforms to an asymmetric Doniach-Sunjic line shape 
in the metal phase. No contamination feature at higher 
binding energy to the main peak is observed. Figure 3(b) 
shows the V 2p core level spectra at 320 and 270 K. The 
V 2p 3 / 2 an d V 2p 1 / 2 main peaks are positioned at ~ 516 
cV and ~ 524 cV. The V 2p spectra show remarkable 
changes across the MIT: the peak shape is sharp in the 
insulating phase, but shows a shoulder structure at low 
binding energy in the metal phase. In addition, the V 
Is core level spectra at 320 and 270 K show very similar 
changes like the V 2p spectra, as shown in Fig. 3(c). 
The V Is main peak is positioned at 5468 cV, with a 
weak satellite at 5481 eV in the insulating phase. The 
main peak and the satellite feature get broadened in the 
metallic phase, with clear additional structure at lower 
binding energies. 

We also carried out cluster calculations using the same 
set of parameters as for the valence band calculations. 
The calculations [Fig. 3(b) and 3(c)] nicely reproduce 
the V 2p and V Is core level experimental spectra. 
The shoulder structure appears in the metallic spectrum, 
whereas the insulating phase spectrum has no shoul- 
der structure in the V Is as well as V 2pi/ 2 and V 
2^3/2 core levels. The results confirm that the screen- 



ing channel from the coherent band at E-p is responsible 
for the shoulder structure. This result is slightly dif- 
ferent from HX-PES results reported in V 2 3 (Refs. H 
and I25T) and Lai-^Sr^MnOs^ which show a clear addi- 
tional peak derived from the coherent screening in both 
HX-PES core level spectra and the cluster model calcu- 
lations. Since the additional peak position is affected by 
the magnitude of the parameters, A, A*, and Udc in the 
calculation^ the additional peak appear well-separated 
from the main peak and the ligand screening (2p 5 3d n L) 
peak in V2O3 and Lai-^Sr^MnOa, resulting in a dis- 
tinguishable peak structure. In contrast, since the A 
is relatively small compared to Udc for VO2, and unlike 
V2O3 and Lai_a;Sr a .Mn0 3 , the coherent screening state 
is positioned near the ligand screening state, resulting 
in a broad shoulder structure. The overall calculation 
results thus reproduce the core level and valence band 
PES spectra and provide consistent electronic structure 
parameters with a U > A. 

In conclusion, we have performed T-dependent valence 
band and core level HX-PES of V02- The valence band 
spectra show gap formation and weight transfer from the 
coherent state (the quasi-particle band, d 2 C — > d 1 C char- 
acter) to a mixed incoherent state [the lower Hubbard 
band, (d 1 — ► d°) + (d 2 L — ► d 1 L) character], supporting a 
MH transition picture. Resonant-PES across the V 2p- 
3d threshold finds that the d}-L state, which overlaps the 
O 2p band, resonates as well as d° and d 1 Q. The spec- 
tral shape changes in V Is and V 2p core levels as well 
as the valence band are nicely reproduced from a cluster 
model, providing a consistent set of electronic structure 
parameters. 

The HX-PES experiments reported here have benefited 
tremendously from the efforts of Dr. D. Miwa of the co- 
herent x-ray optics laboratory RIKEN/SPring-8, Japan 
and we dedicate this work to him. 
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